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Introduction 

Many problems in the optimization theory and other fields of mathe- 
matics are nonlinear in the traditional sense but appear to be linear 
over semirings with idempotent addition. 1 This approach is developed 
systematically as Idempotent Analysis or Idempotent Mathematics (see, 
e.g., [1]— [8]). In this paper we present an idempotent version of Interval 
Analysis (its classical version is presented, e.g., in [9]-[12]) and discuss 
applications of the idempotent matrix algebra to discrete optimization 
theory. 

The idempotent interval analysis appears to be best suited for treat- 
ing problems with order-preserving transformations of input data. It 
gives exact interval solutions to optimization problems with interval un- 
certainties without any conditions of smallness on uncertainty intervals. 
Solution of an interval system is typically ./VP-hard in the traditional 
interval linear algebra; in the idempotent case it is polynomial. The 

* The work was supported by the RFBR grant No. 99-01-01198 and the Erwin 
Schrodinger Institute for Mathematical Physics. 
Submitted to Reliable Computing. 

1 One of the most important examples of an idempotent semiring is the set 
Rmax — K U { — oo} with operations © = max, = + (see example 1.1 in 
subsection 1.2), sometimes called the max-plus algebra. 
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idempotent interval analysis is particularly effective in problems that 
are strongly nonlinear in the traditional sense but possess a hidden 
linear structure over a suitable idempotent semiring, which is often the 
case in optimization theory (see examples in subsection 2.2 below). 

In general, there exists a correspondence between interesting, useful, 
and important constructions and results concerning the field of real 
(or complex) numbers and similar constructions dealing with various 
idempotent semirings. This correspondence can be formulated in the 
spirit of the well-known N. Bohr's correspondence principle in Quantum 
Mechanics; in fact, the two principles are intimately connected (see 
[4, 5, 6]). In a sense, the traditional Mathematics over numerical fields 
can be treated as a quantum theory [4, 5, 13, 14], while the Idempotent 
Mathematics can be treated as a 'classical shadow (or counterpart)' of 
the traditional one. 

In Quantum Mechanics the superposition principle means that the 
Schrodinger equation (which is basic for the theory) is linear. Similarly 
in Idempotent Mathematics the (idempotent) superposition principle 
means that some important and basic problems and equations (e.g., 
the Hamilton-Jacobi equation, which is basic for Classical Mechanics, 
optimization problems, the Bellman equation and its versions and gen- 
eralizations), which are nonlinear in the usual sense, can be treated as 
linear over appropriate idempotent semirings, see [l]-[4]. 

Note that numerical algorithms for infinite-dimensional linear prob- 
lems over idempotent semirings (e.g., idempotent integration, integral 
operators and transformations, the Hamilton-Jacobi and generalized 
Bellman equations) deal with the corresponding finite-dimensional ap- 
proximations. Thus idempotent linear algebra is the basis of the idem- 
potent numerical analysis and, in particular, the discrete optimization 
theory. 

B. A. Carre [15, 16] (see also [17]) used the idempotent linear algebra 
to show that different optimization problems for finite graphs can be 
formulated in a unified manner and reduced to solving Bellman equa- 
tions, i.e., systems of linear algebraic equations over idempotent semir- 
ings. He also generalized principal algorithms of computational linear 
algebra to the idempotent case and showed that some of these coincide 
with algorithms independently developed for solution of optimization 
problems [15]; for example, Bellman's method of solving the shortest 
path problem corresponds to a version of Jacobi's method for solving 
a system of linear equations, whereas Ford's algorithm corresponds to 
a version of Gauss-Seidel's method. 

We stress that these well-known results can be interpreted as a 
manifestation of the idempotent superposition principle. 
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Idempotent Mathematics appears to be remarkably simpler than its 
traditional analog. For example, in the traditional interval arithmetic 
multiplication of intervals is not distributive with respect to addition of 
intervals, while in idempotent interval arithmetic this distributivity is 
conserved. Moreover, in the traditional Interval Analysis the set of all 
square interval matrices of a given order does not form even a semigroup 
with respect to matrix multiplication: this operation is not associative 
since distributivity is lost in the traditional interval arithmetic. On the 
contrary, in the idempotent case associativity is conserved. Finally, in 
the traditional Interval Analysis some problems of linear algebra, such 
as solution of a linear system of interval equations, can be very difficult 
(generally speaking, they are ./VP-hard, see [18]-[24] and references 
therein). We shall show below that in the idempotent case solving an 
interval linear system requires a polynomial number of operations (sim- 
ilarly to the usual Gauss elimination algorithm). Two properties that 
make the idempotent interval arithmetic so simple are monotonicity of 
arithmetic operations and positivity of all elements of an idempotent 
semiring. 

A heuristic introduction into Idempotent Mathematics and a dis- 
cussion of its relations to scientific computing is presented in [4, 13, 
14]); the present paper is concentrated on more technical aspects of 
Idempotent Mathematics, in particular on idempotent matrix theory 
and idempotent interval arithemtics, and their applications to discrete 
optimization. 

This paper consists of five sections. 

The first section concerns general concepts of Idempotent Mathe- 
matics. Subsection 1.1 contains definitions of basic concepts of idem- 
potent arithmetic. Several important examples are presented in sub- 
section 1.2. 

In the second section we present some material from the idempotent 
matrix theory and show its relations to discrete optimization problems. 

Third and fourth sections are central in this paper. The third section 
is devoted to the idempotent interval analysis. To construct an analog of 
calculus of intervals in the context of Idempotent Analysis, we develop 
a set- valued extension of idempotent arithmetic (see subsection 3.1). 
Interval extensions of idempotent semirings are constructed in subsec- 
tions 3.2 and 3.3. In subsection 3.4 we present an idempotent analog 
of E. Kaucher's generalized interval arithmetic. 

In the fourth section we apply the above theory to the problem of 
solving the interval discrete stationary Bellman equation. In particular, 
we discuss construction of an exact outer estimate of a solution set in 
polynomial time and a spectral criterion of convergence of an iterative 
method of solution. 
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In the fifth section we outline a generalization of the idempotent 
interval analysis to the case of positive semirings. 

Some of the results presented in this paper were announced earlier 
in [25, 26]. 

1. Idempotent semirings 

1.1. Basic definitions 

Consider a set S equipped with two algebraic operations: addition © 
and multiplication 0. The triple {S, ©,0} is a semiring if it satisfies 
the following conditions (here and below, the symbol * denotes any of 
the two operations ©, ©): 

• the addition © and the multiplication are associative: x*(y*z) = 
(x*y) * z for all x,y, z G S; 

• the addition © is commutative: x © y = y © x for all x, y G S; 

• the multiplication is distributive with respect to the addition ©: 

x (y © z) = (x y) © (x z) and (x © y) z = (x z) © (y z) 
for all x,y,z G S. 

A semiring S is called idempotent if x © x = x for all x G S. In the 
rest of this paper we shall sometimes drop the word 'idempotent' when 
the corresponding context is clear. 

A unity of a semiring S is an element 1 G S such that 

10X=X01=X 

for all x G S. 

A zero of a semiring S is an element G S such that 0^1 and 
0©x = x, O0x = x0O = O 

for all x G S. 

It is readily seen that if a semiring S contains a unity (a zero), then 
this unity (zero) is determined uniquely. 

A semiring S is said to be commutative if x y = y x for all 
x,yeS. 

Note that different versions of this axiomatics are used; see, e.g., 
[2]-[8] and some literature indicated in these books and papers. 

The addition © defines a canonical partial order on an idempotent 
semiring S: by definition, x =4 y iff xffiy = y. We use the notation x -< y 
if x =<! y and x ^ y. If S contains zero 0, then is its least element 
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with respect to the order =<I. The operations © and are consistent 
with the order ^ in the following sense: if x ==! y, then x*z =^ y*z and 
z-k x =^ z-ky for all x, y, z G 5. 

An idempotent semiring S is said to be a-complete if for any subset 
{x a } C S, including 0, a sum 0{x Q } = Q x a is defined in such a 
way that (0 a x a ) y = Q (x Q y) and y (0 a x a ) = a (y x a ) 
for all y £ S. An idempotent semiring S containing zero is said 
to be b-complete if the conditions of a-completeness are satisfied for 
any nonempty subset {x a } C S that is bounded from above. Any b- 
complete semiring either is a-complete or becomes a-complete if the 
greatest element oo = sup S is added; see [5, 6] for details. 

Note that a x a = sup{x Q } with respect to the canonical partial 
order =<!; in particular, an a-complete idempotent semiring always con- 
tains zero = 00. In an a-complete (6-complete) semiring S the 
inequality a x a ^ a y a holds for all (bounded from above if S is b- 
complete) sets {x a } and {y a } parametrized in such a way that x a ^ y a 
for all a. 

A semiring S with zero is entire if x y = implies that either 
x = or y = for all x, y G S [27]. A semiring 5 is said to satisfy 
the cancellation condition if for all x, y, z G S the equality y = z holds 
whenever x is nonzero and xQy = xQzoryQx = zQx. If a semiring 
satisfies the cancellation condition, then it is entire. A commutative 
semiring S is said to be a semifteld if every nonzero element of S is 
invertible; in this case the cancellation condition holds. 

A semiring S is said to be algebraically closed if the equation x n = y, 
where x n = x • • • x (n times), has a solution for all y G S and 
n G N [28, 29]. Note that in [29] the property of algebraic closedness 
was incorrectly called 'algebraic completeness' due to a translator's 
mistake. 

An idempotent semiring S with zero and unity 1 satisfies the sta- 
bilization condition if the sequence x n © y stabilizes whenever x ^ 1 
and y^O (i.e., x n © y = x n ° © y if n ^ uq for some no) [28, 29]. 

Remark 1.1. In many idempotent semirings algebraic computa- 
tions are greatly simplified by an equality 

(x © y) n = x n © y n . 

For instance, this equality holds in all commutative idempotent semir- 
ings satisfying the cancellation condition (see, e.g., [29], assertion 2.1) 
and in particular in any idempotent semifield. 
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1.2. Examples of idempotent semirings 

The following three examples of idempotent semirings are the most 
important in Idempotent Mathematics. 

Example 1.1. Denote by R max the set S = R U {-oo} equipped 
with operations © = max and = +, where = — oo, 1 = 0. 

Example 1.2. Denote by R min the set R U {+00} equipped with 
operations © = min and = +, where = +00 and 1 = 0. 

Example 1.3. Consider also the set RU {—00, +00} with the oper- 
ations © = max and © = min, where = —00, 1 = +00. 

We see that R max is a 6-complete algebraically closed idempotent 
semifield satisfying the stabilization condition. The idempotent semir- 
ing R m i n is isomorphic to R max . Note that both R max and R m in are 
linearly ordered with respect to the corresponding addition operations; 
the canonical order =<! in R max coincides with the usual linear order ^ 
in R and is opposite to the canonical order =<! in R m i n . 

Example 1.4. Consider the set R max = R max U{oo} with operations 
©, extended by oo©x = 00 for all x € R m ax, oo©x = 00 if £ / and 
00 = 0. It is easily shown that this setjs an a-complete idempotent 
semiring and 00 is its greatest element (R m a X is not a semifield since 
00 is not invertible). 

Example 1.5. Note that the Boolean algebra Sb = {0,1} is a 
unique a-complete idempotent semifield. 

We stress that the equality (x©y) n = x n ®y n holds in the semirings 
of all above examples even though semirings of examples 1.3 and 1.4 
do not satisfy the cancellation condition. 

Example 1.6. Suppose S is an idempotent semiring and X is an ar- 
bitrary set. The set Map(X; S) of all functions X — > S is an idempotent 
semiring with respect to the following operations: 

{f ®9)(x) = f(x)®g(x), (/0y)(x) = /(x)0y(x), x e X. 

If S contains zero and/or unity 1, then the functions o(x) = for all 
x 6 X, e(x) = 1 for all x € X are zero and unity of the idempotent 
semiring Map(X; S). It is also possible to consider various subsemirings 
of Map(X;S). 

Let {S*i, S*2, . . . } be a collection of (idempotent) semirings. There 
are several ways to construct a new idempotent semiring derived from 
the semirings of this collection. 

Example 1.7. Let Si be entire idempotent semirings with opera- 
tions ©j, Oj and zeros 0j, i = 1, . . . , n. The set S = (S± \ {0i}) x 
• • • x (S n \ {0 n }) U {0} is an idempotent semiring with respect to the 
following coordinate-wise operations: 

x*y = (x 1 ,...,x n )-k(y 1 ,...,y n ) = (xi *i yi, . . . , x n * n y n ); 
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the element is zero of this semiring. 

Example 1.8. Note that the direct product Si x • • • x S n is also 
an idempotent semiring with respect to the coordinate-wise opera- 
tions, even if primitive semirings are not entire; its zero is the element 
(Oi,... ,0 n ). 

Note also that even if primitive semirings Si in examples 1.6—1.8 
are linearly ordered sets with respect to the orders induced by the 
correspondent addition operations, the derived semirings are only par- 
tially ordered. On the other hand, if in examples 1.6-1.8 the equality 
(x © y) n = x n © y n holds in all primitive semirings S, Si, then it holds 
in the derived semirings of these examples as well since the operations 
in these semirings are pointwise. 

Many additional examples can be found, e.g., in [2]-[8]. 



2. Idempotent matrices and optimization on graphs 

2.1. Generalities 
2.1.1. Basic definitions. 

Let S be an idempotent semiring. Denote by M&t mn (S) the set of 
all matrices with m rows and n columns whose coefficients lie in an 
idempotent semiring S. 

The sum © of matrices A = (a^), B = (pij) G M.at mn (S) can be 
defined as usual: 

A © B = (aij © bij) G Mat mn (S). 

Let be the corresponding canonical order on the set Mat mn (<S r ). 

The product of two matrices A G Matj m (S) and B G Mat mn (S l ) is 
the matrix 

AB = ^0 a ik b k ^j G Mat fa (5). 

LEMMA 1. The matrix multiplication is consistent with the canon- 
ical order =<! in the following sense: for all A = (aik), A' = ( a 'ik) ^ 
Mat lm (S), B = (b kj ),B' = (b' kj ) G Mat mn (5) ; if A 4 A' in Mat, m (S) 
and B ^ B' in M&t mn (S), then AB A'B' in Mat fa (S'). 
Proof. 

(m \ / m \ / m \ 

a lk b kj U a ik b' kj U a' ik b' kj = A'B', 
k=i / \k=i / \k=i / 

since the operations © and are consistent with the canonical order 
in S. 
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It is easily checked that the set Mat nn (S) of square matrices of order 
n is an (in general, non-commutative) idempotent semiring with respect 
to these operations. Note that we make a slight abuse of notation when 
denote multiplication in this semiring by AB instead of AQ B. 

If is zero of S, then the matrix O = (o^), where Oij = 0, i,j = 
1, . . . , n, is zero of Mat nn (S'); if 1 is unity of S, then the matrix E = 
(5ij), where 5ij = 1 if i = j and 5ij = otherwise, is unity of Mat nn (S). 

A straightforward calculation shows also that if a semiring S is a- 
complete (respectively, 6-complete), then Mat nn (5) is an a-complete 
(respectively, 6-complete) semiring for all n ^ 1. On the contrary, the 
multiplication operation in a matrix semiring Mat nn (S l ) is noncommu- 
tative and does not satisfy the cancellation condition even if the scalar 
multiplication in the semiring S has these properties. 

Let us remember that in the traditional mathematics matrices are 
a kind of coordinate notation for linear operators acting in finite- 
dimensional linear spaces. An obvious analog of the notion of linear 
space in Idempotent Mathematics is the notion of semimodule over 
an idempotent semiring or semifield. In particular, the direct product 
S x • • • x S = S n with coordinate-wise operations of addition and 
multiplication by a scalar from S can be considered a finite-dimensional 
idempotent linear space. Now the correspondence between matrices 
from Mat mn (S l ) and linear operators (homomorphisms) acting from S n 
to S m can be established in the standard way. In particular, Mat nn (<S') 
corresponds to the semiring of endomorphisms of S n . 

Note also that M&t mn (S) itself becomes a linear space if the multi- 
plication by a scalar c G S is defined by c A = (c0Oij) G M&t mn (S) 
for each A = (a^) G M&t mn (S). 



2.1.2. Matrices and graphs. 

Suppose that S is a semiring with zero and unity 1. It is well-known 
that any square matrix A = (a^) G Mat nn (S) specifies a weighted 
directed graph. This geometrical construction includes three kinds of 
objects: the set X of n elements xi,...,x n called nodes, the set T of all 
ordered pairs (xi,Xj) such that a^- ^ called arcs, and the mapping 
A: r — > S such that A(xi,Xj) = aij. The elements a^- of the semiring 
S are called weights of the arcs. 

Conversely, any given weighted directed graph with n nodes specifies 
a unique matrix A G M&t nn (S). 

This definition allows for some pairs of nodes to be disconnected if 
the corresponding element of the matrix A is and for some channels 
to be 'loops' with coincident ends if the matrix A has nonzero diagonal 
elements. This concept is convenient for analysis of parallel and dis- 
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tributed computations and design of computing media and networks 
(see, e.g., [30, 31]). 

Recall that a sequence of nodes of the form 

P = (yo,Vi, ■■■ ,Vk) 

with k ^ and (yi, yi+i) G I\ i = 0, . . . , k — 1, is called a path of length 
k connecting yo with y k . Denote the set of all such paths by Pk(yo, yu)- 
The weight A{p) of a path p € Pk(yo,yk) is defined to be the product 
of weights of arcs connecting consecutive nodes of the path: 

A(p) =A(y ,y 1 )Q---QA(y k - 1 ,y k ). 

By definition, for a 'path' p £ Po(xi,Xj) of length k = the weight is 
1 if i = j and otherwise. 

For each matrix A <G Mat nn (5) define A = E = (<5y) (where 5ij = 1 

if % = j and <% = otherwise) and A k = AA k ~~ l , k ^ 1. Let a-^ be the 
(i, j)th element of the matrix A fc . It is easily checked that 

Oj n © ' ' ' © a ik-iik ■ 

io=i, ik=3 

Thus is the supremum of the set of weights corresponding to all 
paths of length k connecting the node Xi = Xi with Xi k = Xj. 

2.2. Matrix formulation of some optimization problems 

2.2.1. Closure operation and the algebraic path problem 
Suppose S is an idempotent semiring with unity 1. The closure opera- 
tion * in S is defined by a 'power series' 

x* = i e x e x 2 e • • • 

for any x £ S. This operation was first introduced by S. Kleene in a 
special case [32]; it is well-known in the context of Idempotent Analysis 
[7, 8, 15, 16, 17]. 

Of course, the sum of this power series must be well-defined. In 
particular, the infinite sum (Bo<fc<oo %k ^ s defined in every a-complete 
semiring. In the semirings of examples 1.1 and 1.2 the closure x* is de- 
fined for all x such that x ^ 1 (so x* = 1). In the (a-complete) semirings 
of examples 1.3-1.5 the closure is defined for all their elements. 

LEMMA 2. The closure operation is consistent with the canonical or- 
der =4 in S in the following sense: if x,x' G S and x =4 x' , then 

x* 4 {x'y. 



M _ 
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Proof. Since the operation is consistent with the canonical order 
in S, the inequality x k ^ (x') k holds for all k ^ 0. Thus x* = 

In the matrix semiring Mat nn (5) the closure is defined by 

A* = E@ A® A 2 . 

(*) 

Denote the elements of the matrix A* by a- , i, j = 1, . . . , n; then 

«!? = © © A <p)- 

0^fc<oo p&Pk(xi,Xj) 

The closure matrix ^4* solves the well-known algebraic path prob- 
lem, which is formulated as follows: for each pair (xi,xj) calculate the 
supremum of weights of all paths (of arbitrary length) connecting node 
Xi with node Xj. The closure operation in matrix semirings has been 
studied extensively (see, e.g., [3]-[17] and references therein). 

Example 2.1. The shortest path problem. Let S = IW n , so 
the weights are real numbers. In this case 

A(p) = A(y ,y 1 ) + A(y 1 ,y 2 ) H h A(y k - 1 ,y k ). 

If the element cijj specifies the length of the arc (x^ Xj) in some metric, 

<*) 

then a\j is the length of the shortest path connecting x i with Xj . 

Example 2.2. The maximal path width problem. Let S = 
M. U {0, 1} with © = max, = min as in example 1.3. Then 

a a = , , ma ? ,^(p)> a (p) = min (^(yo,yi), •• ■ ,Mvk-i,yk))- 

pe U Pk{xi,xj) 

If the element specifies the 'width' of the arc (xj, Xj), then the width 
of a path p is defined as the minimal width of its constituting arcs and 
the element a\j gives the supremum of possible widths of all paths 
connecting X{ with Xj. 

Example 2.3. A simple dynamic programming problem. Let 
S = M max and suppose aij gives the profit corresponding to the tran- 
sition from Xi to xj. Define the vector B = (pi) 6 Mat„i(R max ) whose 
element bi gives the terminal profit corresponding to exiting from the 
graph through the node X{. Of course, the negative profits (or, rather, 
losses) are allowed. Let m be the total profit corresponding to a path 
p G P k (xi,Xj), i.e. 

m = A(p) + bj. 

Then it is easy to check that the supremum of profits that can be 
achieved on paths of length k beginning at the node x% is equal to 
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{A k B)i and the supremum of profits achievable without a restriction 
on the length of a path equals {A*B)i. 

Example 2.4. The matrix inversion problem. Note that in the 
formulas of this section we are using distributivity of the multiplication 
with respect to the addition © but do not use the idempotency axiom. 
Thus the algebraic path problem can be posed for a nonidempotent 
semiring S as well (see, e.g., [33]). For instance, if S = M, then 

A* = E + A + A 2 + • • • = (E - A)- 1 . 

If \\A\\ > 1 but the matrix E — A is invertible, then this expression 
defines a regularized sum of the divergent matrix power series X^«>o ^ ' 

2.2.2. Discrete stationary Bellman equation. 

Bellman, Isaacs, and Hamilton-Jacobi equations are central in different 
parts of optimization theory. It is well-known that these equations are 
strongly nonlinear in the traditional sense but have a linear structure 
over appropriate idempotent semirings (see, e.g., [l]-[7]). This fact 
was first observed by B. A. Carre for discrete versions of Bellman 
equation [15]-[17]. 

The following equation (the discrete stationary Bellman equation) 
plays an important role in both discrete optimization theory and idem- 
potent matrix theory: 

X = AX © B, 

where A £ Mat nn (S), X, B € Mat ns (S l ); matrices A, B are given and 
X is unknown. The discrete stationary Bellman equation is a natural 
counterpart of the usual linear system AX = B in traditional linear 
algebra. 

Note that if the closure matrix A* = E © A © A 2 © • • • exists, 
then the matrix X = A*B satisfies the discrete stationary Bellman 
equation because A* = AA* © E. It can be easily checked that this 
special solution is the minimal element of the set of all solutions to the 
discrete stationary Bellman equation. 

We emphasize that this connection between the matrix closure op- 
eration and solution to the Bellman equation gives rise to a number 
of different algorithms for numerical calculation of the closure matrix. 
All these algorithms are adaptations of the well-known algorithms of 
the traditional computational linear algebra, such as Gauss-Jordan 
elimination, various iterative and escalator schemes, etc. 

In fact, the theory of the discrete stationary Bellman equation can 
be developed using the identity A* = AA* © £ as an additional axiom 
without any substantive interpretation (the so-called closed semirings; 
see, e.g., [17, 34]). 
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2.3. TWO KNOWN THEOREMS 

In this subsection we recall some general results of the idempotent 
matrix theory that are necessary for the subsequent sections. 

2.3.1. Existence of a closure 

Suppose A £ Mat nn (S') and an idempotent semiring S is not a-complete. 
Then a closure matrix A* = E © A © A 2 © • • • might not be defined 
if this series diverges. Let us formulate a sufficient condition for the 
existence of a closure, following the work of B. A. Carre [15]. 

A matrix A = (a^) € M&t nn (S) is said to be definite (respectively, 
semi- definite) if 

A(p) -< 1 (respectively, A(p) =^ 1) 

for any path p £ Pk(yo, Uk) such that yo = yk, k ^ 1 (i.e., for any closed 
path). Obviously, every definite matrix is semi-definite. This definition 
is similar to that of [15] but B. A. Carre considers an ordering that is 
opposite to =<;. 

THEOREM 1. (Carre). Let A be a semi-definite square matrix of order 
n. Then 

k n—l 
1=0 1=0 

for n—l, so the closure matrix A* = ©^ ^4 fe exists and is equal 
to®t = lAK 

For the proof see, e.g., [15], Theorem 4.1. The basic idea of the proof 
is evident: in the graph of a semi-definite matrix it is impossible to 
construct a path of arbitrarily large weight since the weight of any 
closed part of a path cannot be greater than 1. Thus there exists a 
universal bound on path weights, which makes truncation of the infinite 
series expressing the closure matrix possible. 

2.3.2. Eigenvectors and eigenvalues 

The spectral theory of matrices whose elements lie in an idempotent 
semiring is similar to the well-known Perron-Frobenius theory of non- 
negative matrices (see, e.g., [3, 8, 28, 29]). 

Recall that a matrix A = iaij) £ Mat nn (S') is said to be irreducible 
in the sense of [8] if for any 1 ^ i,j ^. n there exist an integer k ^ 1 
and a path p G Pk(xi, Xj) such that A(p) ^ 0. In [28, 29] matrices with 
this property are called indecomposable. 

We borrow the following important result from [28, 29] (see also [8]): 
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THEOREM 2. (Dudnikov, Samborskii). // a commutative idempotent 
semiring S with a zero and a unity 1 is algebraically closed and sat- 
isfies the cancellation and stabilization conditions, then for any matrix 
A G M&t nn (S) there exist a nonzero 'eigenvector' V G Mat n i(5) and 
an 'eigenvalue' A G S such that AV = A V. If the matrix A is 
irreducible, then the 'eigenvalue ' A is determined uniquely. 

For the proof see [29], Theorem 6.2. 

3. Idempotent interval arithmetics 

3.1. Set-valued idempotent arithmetics 

Suppose S is an idempotent semiring and S is a system of its subsets. 
Denote the elements of S by x, y, . . . Recall that the symbol * denotes 
any of the operations ©, in the semiring S (see section 1.1). Define 
x*y = | x G x,y G y }. 

We shall suppose that S satisfies the following two conditions: 

1. If x, y G S and * is an algebraic operation in S, then there exists 
z G S such that z D x * y. 

2. If {z a } is a subset of S such that f] z a ^ 0, then there exists the 
infimum of {z a } in S with respect to the ordering C, i.e., the set 
x G S such that x C f] a z a and y C x for any y G S such that 

Define algebraic operations ©, in 5 as follows: if x, y G S, then 
x*y is the infimum of the set of all elements z G S such that z D x*y. 
Thus x*y is 'the best upper estimate' for the set x*y in S. 

PROPOSITION 1. The following assertions are true: 

• S is closed with respect to the operations ©, 0. 

• // the system S contains all one-element subsets of S, then the 
semiring {S, ©, 0} is isomorphic to a subset of the algebra {S, ©, ©}. 

The proof is straightforward. 

The following example shows that not much can be said in general 
about the algebra {<S, ©, 0}. 

Example 3.1. Let S = 2 s ; thus x*y G S for all x, y G S, 
so x ★ y = x * y. In general, the set S with these 'naive' operations 
©, satisfies the above assumptions but is not an idempotent semiring. 
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Indeed, let S be the semiring (R max \ {0}) x (R max \ {0}) U {0} with 
coordinate- wise operations ffi, (see example 1.7). Consider a set x = 
{(0, 1), (1, 0)} G S; we see that 

xex = {(o,i),(i,o),(i,i)}/x 

and if y = {(1, 0)}, z = {(0, 1)}, then 

x©(y©z) = {(1, 2), (2, 1)} + (x©y)©(x0z) = {(1, 1), (1, 2), (2, 1), (2, 2)} 

This means that the system S with the operations ©, does not satisfy 
axioms of idempotency and distributivity. 

It follows that S should satisfy some additional conditions in order 
to have the structure of an idempotent semiring. In the next sections 
we consider the case when S is a set of all closed intervals; this case 
is of particular importance since it represents an idempotent analog of 
the traditional Interval Analysis. 

3.2. Weak interval extensions of idempotent semirings 

Let S be a set partially ordered by a relation =<!. A (closed) interval in 
S is a subset of the form x = [x,x] = {t e 5 | t ^ x}, where x, 
x G S (x =^ x) are called the lower and the upper bound of the interval 
x, respectively. 

Note that if x and y are intervals in S, then xCyiffy^x^x^y. 
In particular, x = y iff x = y and x = y . 

Example 3.2. Let x, y be intervals in an idempotent semiring S 
with the canonical partial order =<;. In general, the set x*y is not an 
interval in S. Indeed, consider a set S = {0, a, b, c, d} and let © be 
defined by the following order relation: is the least element, d is the 
greatest element, and a, b, and c are noncomparable with each other. 
If is the zero multiplication, i.e., if x y = for all x, y € S, then S 
is an idempotent semiring without unity. Let x = [0, a] and y = [0, b]; 
then x © y = {0, a, b, d}. This set does not contain c and hence is not 
an interval since ^ c ^ d. 

Let S be an idempotent semiring. We define a weak interval exten- 
sion I(S) of the semiring S to be the set of all closed intervals in S 
equipped with the following operations ffi, 0: x*y = [x*y,x*y] for 
all x,y G I(S), where * denotes or 0. 

PROPOSITION 2. The weak interval extension I(S) of the idempotent 
semiring S is closed with respect to the operations ©, and forms an 
idempotent semiring. 

Proof. The set I(S) with the operations ©, can be identified with a 
subset of an idempotent semiring SxS with coordinate- wise operations 
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(see example 1.8). Since x*y x*y whenever x x and y =<: y , I(S) 
is closed with respect to the operations 0, 0; hence it is an idempotent 
semiring (a subsemiring of S X S) . 

The operation © generates the corresponding canonical partial order 
=^ in I(S): x ^ y iff x < y and x =<! y in S. 

The following proposition shows that this choice of operations 0, 
is consistent with the general construction described in the previous 
section. 

PROPOSITION 3. For all x,y G I(S) the interval x*y contains the 
se£ x*y and is the least interval ofl(S) with this property. In particular, 
bounds of the interval x*y belong to x*y. 

Proof. Let z G I(S) be such that x*y C z. We have x*y G x*y C z 
and x*yGx*yCz; thus z =<! x * y and x ★ y ^ z. This means that 
x ★ y C z, i.e., that the interval x ★ y is contained in any interval 
containing the set x * y. 

Now take t G x * y and let x G x, y G y be such then t = x*y. By 
definition of an interval, x =^ x =^ x and y =<! y =^ y . Since operation ★ 
is consistent with the order ^, we see that x*y^x*y^x*y; this 
means that t G x*y, that isx*yCx*y. This concludes the proof. 

COROLLARY (monotonicity property). I/x C xi, y C yi, thenx*y C 
xi*yi. 

Remark 3.1. Note that in general the system S = I(S) of subsets 
of the semiring S does not satisfy condition 2 of section 3.1 if S is not 
6-complete. 

Let an idempotent semiring S be a-complete (respectively, 6-comple- 
te) and {x Q } be an infinite subset of its weak interval extension I(S) 
(with an additional requirement in the case of 6-complete S that {x a } 
is bounded from above with respect to the canonical order ^ in I(S)). 
Define the (infinite) sum of elements of this subset by 

X « = 

a 

PROPOSITION 4. If the semiring S is a-complete (respectively, b- 
complete), then the semiring I (S) is a-complete (respectively, b-complete) 
with respect to the above definition of an infinite sum. 

Proof. Evidently, the interval Q x Q is well-defined if the subset 
{x a } satisfies the above conditions. Now we shall check the distribu- 
tivity axiom. 



a a 
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If S is a-complete and X C I(S) is empty, then X = [0, 0] and 
y0(ei) = (©*) 0y = [0,0] for all y € I(S). If X = {x a } is 
nonempty and infinite, then by a straightforward calculation 



y© 0x c 



0(yoxj,0(y©x a ) 



0(y©x a ) 



and similarly (0 a x a ) y = a (x a y) for all y € I(S). Thus the 

idempotent semiring I(S) is a-complete (6-complete) if S is a-complete 
(6-complete) . 

In what follows, we shall always assume that all infinite sums in weak 
interval extensions of a-complete and 6-complete idempotent semirings 
are defined as above. 

PROPOSITION 5. The interval Q x Q contains the set Q x Q = 
{ ©a x a \ x a £ x Q for all a } and is the least interval of I(S) with this 
property. In particular, bounds of the interval Q x Q belong to Q x a . 



The proof is similar to the proof of proposition 3. 
The following two propositions are straightforward consequences of 
our definition of the operations ©, 0: 

PROPOSITION 6. // an idempotent semiring S is commutative, then 
the semiring I(S) is commutative. 

PROPOSITION 7. // an idempotent semiring S contains zero (re- 
spectively, unity 1), then the interval [0,0] (respectively, [1,1]) is zero 
(respectively, unity) of I(S). 

PROPOSITION 8. // S is entire an idempotent semiring, then I(S) 
is also entire. 

Proof. Let x, y € I(S) and x / [0, 0], y / [0, 0]. Recall that x =^ x, 
y =<! y ; thus 0, y ^ 0. If z = x y , then z = x0y/O, since S is 
entire. It follows that z / [0,0]. 

PROPOSITION 9. // S is algebraically closed and for all x,y £ S, 
n G N the equality (x © y) n = x n © y n holds, then I(S) is algebraically 
closed. 

Proof. Suppose x ra = x© . . . ©x = y. By definition of the operations 
©, 0, we see that x n = y and x™ = y. Let z G S and z <G S be the 
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solutions of these two equations. We claim that z and z can be chosen 
such that z =<! z, i.e., the interval [z, z] is well defined. 

Take z' = z © z; hence z z'. Since z' n = (z © z) n = z™ z™ in 5, 
we see that z' n = y © y = y . Thus [z, z'} n = [y , y] = y . 

Remark 3.2. Recall that the equality (x © y) n = x n ®y n holds in 
many semirings, including all semirings listed in examples 1.1—1.5. 

3.3. Strong interval extension 

We stress that in general a weak interval extension I(S) of an idem- 
potent semiring 5 with zero and unity 1 that satisfies cancellation 
and stabilization conditions does not inherit the latter two properties. 
Indeed, let x z = y z, where z = [0,z] and z / 0; then z is a 
nonzero element but this does not imply that x = y since x and y may 
not equal each other. Further, let y = [0,y] ^ [0,0]; then the lower 
bound of x™0y may not stabilize when n — > oo. 

Therefore we define a strong interval extension of an idempotent 
semiring 5 with zero to be the set 1(5) = {x = [x,x] € 7(5) | -< 
x =<! x } U {[0, 0]} equipped with operations ©, defined as above. It 
is clear that 1(5) C I(S). 

Note that this object may not be well-defined if the semiring 5 is 
not entire. Indeed, let intervals x, y G 1(5) be such that -< x -< x, 
-< y -< y, x©y = 0, and x0y ^ 0; then x0y = [O,x0y] £ 1(5). 

Throughout this section, we will suppose that the strong interval 
extension 1(5) of an idempotent semiring 5 is closed with respect to 
the operations © and 0. To achieve this, it is sufficient to require that 
the semiring 5 is entire. 

THEOREM 3. The strong interval extension 1(5) of an idempotent 
semiring S is an idempotent semiring with respect to the operations © 
and with zero [0,0]. It inherits some special properties of the semiring 
S: 

• If S is a-complete (respectively, b-complete), thenl(S) is a-complete 
(respectively, b-complete). 

• If S is commutative, then 1(5) is commutative. 

• If 1 is unity of S, then [1, 1] is unity o/I(5). 

• If S is entire, then 1(5) is entire. 

• If S is algebraically closed and for all x,y € S, n G N the equality 
[x © y) n = x n © y n holds, then 1(5) is algebraically closed. 
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• If S satisfies the cancellation condition, then 1(5) satisfies the 
cancellation condition. 

• If S is a semiring with unity 1 satisfying the stabilization condi- 
tion, then the semiring 1(5) satisfies the stabilization condition. 

Proof. Using the definition of the operations ©, and proposition 2, 
it is easy to check that 1(5) is an idempotent semiring with respect to 
the operations ©, © and contains zero element [0,0]. Propositions 4 
and 6-9 imply the first five statements. 

Suppose 5 satisfies the cancellation condition, x, y, z € 1(5), and z 
is nonzero. Ifx©z = y©z, then x z = y z and x z = y z; since 
z ^ [0, 0] in 1(5), and z / 0, and it follows from the assumptions 
that x = [x, x] = [y, y] = y. If z x = z y , then x = y similarly. 

Suppose further that 5 satisfies the stabilization condition; by defini- 
tion, y ^ and y / for any nonzero y £ 1(5). Consider the sequence 
x n ©y; stabilization holds in 5 for both bounds of the involved intervals 
and hence, by definition of the operations ffi, 0, for the whole intervals 
as elements of 1(5). 

Suppose 5 is an idempotent semiring; then the map t: 5 — > I(S) 
defined by i{x) = [x, x] for all x £ 5 is an isomorphic imbedding of 5 
into its weak interval extension I(S). If the semiring 5 has zero and 
its strong interval extension 1(5) is well-defined, then the map t takes 
5 in 1(5) C I(S), so it is an isomorphic imbedding of 5 into its strong 
interval extension. To simplify notation in the sequel, we will identify 
the semiring 5 with subsemirings l(S) C I(S) or l(S) C 1(5) C I(S) 
and denote the operations in 7(5) or 1(5) by ©, 0. If the semiring 5 
contains unity 1, then we denote the unit element [1, 1] of I(S) or 1(5) 
by 1; similarly, we denote [0,0] by 0. 

Also, we shall drop the word 'strong' and call 1(5) simply an 'interval 
extension' of the semiring 5. 

3.4. An idempotent analog of the Kaucher interval 
arithmetic 

We stress that in idempotent interval mathematics most of algebraic 
properties of an idempotent semiring are conserved in its interval ex- 
tension. On the other hand, if 5 is an idempotent semifield, then the set 
1(5) is not a semifield but only a semiring satisfying the cancellation 
condition. 

Recall that any commutative idempotent semiring 5 with a zero 
can be isomorphically embedded into an idempotent semifield 5 
provided that 5 satisfies the cancellation condition (see, e.g., [28]). 



Rc2_l.tex; 1/02/2008; 6:57; p. 18 



Idcmpotent Interval Analysis and Optimization Problems 



19 



If S coincides with its subsemifield generated by S, then S is called 
a semifield of fractions corresponding to the semiring S. Consider the 
following equivalence relation: for any (x,y), (z,t) <G S x (S \ {0}) 

(x,y) ~ (z,t) iff xQt = yOz. 

Then the semifield of fractions can be constructed as the quotient S x 
(S\ {0})/ ~, equipped with operations 

(x,y)®(z,t) = ({xQt)®{yQz),yQt), (x,y)Q(z,t) = (xQz,yQt). 

The pairs (x, y) behave as 'fractions' with the 'numerator' x and the 
(nonzero) 'denominator' y with respect to the above operations. It is 
easy to check that these operations satisfy the axioms of a commutative 
idempotent semiring with a zero element { (0, y) j y ^ } and a unity 
{{y->y) I D 0}. F° r every 'fraction' (x,y) representing a nonzero 
element of S its inverse element is given by the fraction (y,x); hence 
this algebraic structure satisfies all axioms of an idempotent semifield. 

In the context of the traditional Interval Analysis this is close to 
the construction of the Kaucher interval arithmetic [35, 36]. In ad- 
dition to usual intervals [x,y], where x ^ y, this arithmetic includes 
quasi- intervals [x, y] with y ^ x, which arise as inverse elements for 
the former with respect to addition. In contrast, in the idempotent 
case quasi-intervals arise as inverse elements with respect to semiring 
multiplication. 

The following statement shows that in this case the semifield of 
fractions of interval extension 1(5) corresponding to an idempotent 
semiring S with cancellation condition has a very simple structure: it 
is isomorphic to the idempotent semifield (S'\{0}) x (S'\{0})U{(0, 0)} = 
(S\ {0}) 2 U {(0,0)} (see example 1.7). 

PROPOSITION 10. Suppose S is a commutative idempotent semir- 
ing with a zero 0, S satisfies the cancellation condition, and S is its 
semifield of fractions; then a semifield of fractions corresponding to 
the interval extension I(S) is isomorphic to the semifield (S \ {0}) 2 U 
{(0,0)} with coordinate-wise operations. 

Proof. It follows from theorem 3 that I(S') is a commutative idem- 
potent semiring with a zero element = [0, 0] and satisfies the can- 
cellation condition. Thus 1(S) can be isomorphically imbedded into its 
semifield of fractions. 

Define the map ip: I(S) x (I(S) \ {0}) -» (S \ {0}) 2 U {(0,0)} by 
the rule </?((x,y)) = (x y^ 1 ,^ y" 1 ), where inverses are taken in 

the semifield S. This map is surjective. Indeed, (0,0) = </?((0,y)) for 
any y / 0; let us check that if a, b G S, a / 0, b ^ 0, then there 
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exist x, y G y / 0, such that (a,b) = y((x, y)). By definition of 

a semiring of fractions, there exist nonzero a±, a 2 , b±, b 2 £ 5 such that 
a = ai a^ 1 , b = b\Q b^ 1 in S. Define 

x = ai 61 b 2 , x = (ai 61 6 2 ) © (a 2 6?), 
y = a 2 61 6 2 , y = (ai b 2 2 ) (a 2 61 b 2 ); 

thus O^x^x, 0-<y=^y and 

x y" 1 = 01 61 6 2 &2 1 © © a 2 1 = a i © a 2 1 = a ' 
x y _1 = 61 (ai 6 2 © a 2 61) (ai b 2 © a 2 63 1 = b - 

Since x0y _1 = z0t _1 iff xQt = y Qz for any x,y, z,t € 5 such that 
1/ / and i / 0, we see that <p((x, y)) = </?((z,t)) iff (x, y) ~ (z,t). 
Also, 

tp((x, y) © (z, t)) = p(((x t) © (y z), y t)) 

= ((x y" 1 )) © (z t" 1 ), (x y- 1 )) © (z t" 1 )) 

= v((x,y))©v((z,t)), 
^((x,y)©(z,t)) = ^((x0z,y0t)) 

= ((x y" 1 ) (z t" 1 ), (x y" 1 ) (z r 1 )) 
= ¥>((x,y))0^((z,t)). 

Thus the mapping is an isomorphism of the semifield of fractions 
corresponding to 1(5) and the idempotent semifield (S \ {0}) 2 U {0}. 

The commutativity condition in this proposition is a natural one. 
Indeed, it follows from the theory of ordered groups that if S is a 
6-complete idempotent semiring such that its nonzero elements are 
invertible, then S is commutative and hence is a semifield (see, e.g., 
[5])- 



4. Application to algebraic problems arising in discrete 

optimization 

The discrete stationary Bellman equation and the idempotent matrix 
closure operation are substantial for discrete optimization theory. In 
this section we consider two algebraic problems arising in the case of 
interval Bellman equation: construction of exact interval estimates for 
solution and convergence of an iterative method of solution. 
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4.1. Preliminaries. 



Suppose 5 is an idempotent semiring and I(S) is its weak interval 
extension; then Mat nn (I(S)) is an idempotent semiring. If the interval 
extension 1(5*) of the semiring 5 is well-defined, then Mat nn (I(5)) is 
an idempotent semiring. 

If A = (atj) G Mat mn (/(5)) [A = (a*,) G Mat m „(I(5))] is a (not 
necessarily square) interval matrix, then the matrices A = (a^) and 

A = are called the lower and the upper matrix of the interval 

matrix A. 



PROPOSITION 11. Let S be an idempotent semiring. The mapping 

A £ Mat nn (/(£)) i— ► A, A € 7(Mat„ n (5)) is an isomorphism of 

idempotent semirings Mat nn (/(5)) and I(M&t nn (S)). If the semiring 5 
has an interval extension 1(5), then the mapping A G Mat nn (I(5)) i— > 

A, A G 7(Mat nn (5)) is an isomorphism of idempotent semirings 

Mat nn (I(5)) and I(Mat nn (S)). 



Here intervals 



A, A 



in I(Mat nn (5)) are defined with respect to 

the canonical partial ordering =<! in Mat nn (5) (see section 2.1). The 
proof follows easily from the definition of the operations ©, 0; in- 
deed, this definition implies that addition (respectively, multiplication) 
of interval matrices is reduced to the separate addition (respectively, 
multiplication) of their lower and upper matrices. 

Of course, the notation A G A means that A =<! A =4 A. 

The following proposition is an immediate consequence of theorem 2: 

PROPOSITION 12. If a commutative idempotent semiring S with zero 
and unity 1 is algebraically closed and satisfies cancellation and 
stabilization conditions, then for any matrix A G Mat nn (I(<S)) there 
exist a nonzero 'eigenvector' V G Mat n i(I(S 1 )) and an 'eigenvalue' 
A G 1(5) such that AV = A V. If the matrix A is irreducible, then 
the 'eigenvalue ' A is determined uniquely. 

By definition of the operations ©, 0, AV = A0 V and AV = A0 V. 

4.2. Efficient outer interval estimates for solution sets 
of Bellman equations. 

Consider the following interval discrete stationary Bellman equation 
(see also subsection 2.2.1): 

X = AX © B. 
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Here A 6 Mat nn (I(5)), B £ Mat ns (I(S')), and X is an unknown matrix 
of n rows and s columns. 

Following the discussion in [37], we might consider two different 
notions of a solution to the discrete stationary Bellman equation: 

• The united solution set: 

E(A, B) = { X e Mat ns (5) | X = AX@B for some AeA,5eB). 

• The algebraic solution: X 6 Mat ns (I(5)) such that X = AX © B. 

For definition of some other possible solution sets and discussion of 
their relations see [37, 38]. 

Let us remember that the minimal solution to the equation X = 
AX@B in the sense of the canonical order =<; in Mat ns (5) is X = A*B. 
In what follows, we shall always suppose that the closure matrix exists 
and consider only minimal solutions. Recall that if the matrix A is def- 
inite in the sense of B. A. Carre (see [15] and subsection 2.3), then this 
solution is unique. We shall use the term united minimal solution set 
for the united solution set S(A, B) consisting only of minimal solutions 
and denote it by £ m i n (A, B). 

In the traditional interval analysis the united solution set has a very 
complicated structure and requires exponentially many operations for 
its full description. Even the problems of recognition whether this set 
is empty and finding an outer interval estimate within a given error 
for this set can be iVP-hard (see [18, 19], [20]-[24] and references in 
these papers and surveys; see also a discussion in [37, 38]). However, if 
all interval entries of the interval matrix A consist only of nonnegative 
numbers, the algebraic solution of the system X = AX + B turns to 
be a sharp outer interval estimate of the united solution set [39] (see 
also [11], Theorem 12.2). The following result shows that in the idem- 
potent Interval Analysis the (minimal) algebraic solution A*B of the 
equation X = AX © B is a sharp outer interval estimate of the united 
minimal solution set £ m j n (A, B) for all matrices A € Mat„ n (I(S')), 
B € Mat ns (I(S)). 

THEOREM 4. The interval matrix A*B <= Mat ns (/(£)), considered 
as an element of I(M&t ns (S)), contains the united minimal solution set 
E m j n (A,B) and is the least interval of I(M&t ns (S)) with this property. 
In particular, bounds of the interval A*B belong to S m i n (A,B). 

Proof. Since matrix multiplication and the closure operation are 
consistent with the canonical order ^ in matrix semirings (see subsec- 



tions 2.1.1 and 2.2.1), we see that A*B 
element of I(Mat ns (S)). 



A*B, A B 



is a well-defined 
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Let Z £ Mat ns (I(S')) C I(Mat ns (S)) containthe set £ min (A,B); 
then A*B C Z since A*B £ £ min (A, B) C Z, A*B £ S min (A, B) C Z. 

Conversely, let T £ S min (A, B) and T = A*B, where A £ A, B £ B. 
Then T £ A*B since A*B ^ T A*B. We see that S min (A,B) C 
A*B; this concludes the proof. 

COROLLARY. It is possible to obtain the sharp outer interval estimate 
A*B of the united minimal solution set £ mm (A,B) in a polynomial 
number of operations. 

Proof. By definition of the operations 0, © in the interval extension 
of the idempotent semiring S, operations with interval matrices are 
reduced to separate operations with their lower and upper matrices. 
On the other hand, it is possible to obtain the algebraic solution of 
the discrete stationary Bellman equation X = AX © B by means 
of the Gauss elimination (or some other efficient algorithm of linear 
algebra), which requires a polynomial number of operations. Repeating 
this calculation for lower and upper matrices A*B and A B separately, 
we get the outer interval estimate A*B after a polynomial number of 
operations. 

4.3. Spectral criterion of convergence of iterative 
process. 

Consider the following iterative process: 



where X fc £ Mat ns (I(S)), k = 0, 1, . . . 

The following proposition is due to B. A. Carre [15] (Theorem 6.1) 
up to some terminology: 

PROPOSITION 13. // a matrix A £ Mat„ n (S') is semi-definite, then 
the iterative process X^+i = AX^ ®B stabilizes to the (minimal) solu- 
tion X = A*B of the equation X = AX ffi B after at most n iterations 
for any initial approximation Xq £ Mat ns (S l ) such that X$ A*B. 

Suppose an idempotent semiring S satisfies the assumptions of propo- 
sition 12. Let Ai, . . . , \ q , 1 ^ q ^ n, be the eigenvalues of the matrix 



A £ Mat„ n (I(5)). Denote sup{Ai,... ,A fc } = 0f =1 A; by p(A). It 



is possible to give a simple spectral criterion of convergence of our 
iterative process: 
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THEOREM 5. Let S be a commutative semiring satisfying conditions 
of proposition 12 and matrix A G Mat nn (I(S')) be such that p(A) =^ 1. 
Then the iterative process X^ + i = AX& © B ; k ^ 0, stabilizes to the 
(minimal) algebraic solution X = A*B of equation X = AX©B after 
at most n iterations for any Xo G Mat ns (I(S')) such that Xo ^ X. 

Proof. It follows from the definition of the operations ffi, that it is 
sufficient to prove that sequences of lower and upper matrices of {X^} 
converge separately. To this end, we shall show that the matrices A and 
A are semi-definite; then the result will follow from proposition 13. 

Since a^- =^ aTJ for all we need only to prove that A is semi- 
definite. First we shall prove this if A is irreducible. Using the expres- 
sion for a unique eigenvalue of an irreducible matrix A in terms of cycle 
invariants [29] 

A = [ a ii«2 • • • a iih]^ ^ j 

1=1,... ,n 

(h,-,k) 

where (p(n) is the least common multiple of the numbers 1, . . . , n, we see 
that for any closed path p its weight A(p) = • • • a^ satisfies 
A(p) =<; 1 if A 1 (indeed, if A(jp) ffi 1 ^ 1, then, by remark 1.1, 

(1 © A(p)y( n )l l = le A( P y( n )l l y l, so 1 © A v(n) >- 1, this is a 
contradiction). Thus A is a semi-definite matrix. 

Now let A be a reducible matrix. It follows from idempotent matrix 
algebra (see, e.g., [8]) that there exists a permutation of rows and 
columns of matrix A taking it into an upper block triangular matrix 

(B 1 * 
D _ B 2 



\0 ■•• B k J 

and all square matrices Bi,... ,B k are either zero or irreducible. Every 
eigenvalue A of A is an eigenvalue of B; we claim that in fact it is an 
eigenvalue of some Bi, I = 1, . . . , k. Indeed, let V = (vi) G Mat n i(S) 
be an eigenvector of B with an eigenvalue A. Consider a decomposition 
of the set of nodes X = X\ U • ■ ■ U where X\ n X s = if I ^ s and 
Bi = (bij) XitXjeXl , I = 1, . . . , k; let Z = max{ I \ Vi ^ for some Xi € 
X[ }. We see that A is a unique eigenvalue of the irreducible matrix 
Bi corresponding to an eigenvector (vi) Xie x io of B[. The condition 
p(A) 1 implies that Bi,... ,Bk are semi-definite. Since there is no 
closed path p containing nodes Xj G Xi, Xj G X s , I ^ s, we conclude 
that A is a semi-definite matrix. 



1 < k ^ n, 
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Remark 4.1. Compare this simple proposition with the well-known 
spectral convergence criterion of the iterative process in traditional 
Interval Analysis ([11], Theorem 12.1), which in our notation has the 
following form: 

The iterative process X^ +1 = AX^. + B, k ^ 0, converges to a unique 
solution X of the equation X = AX + B for any Xo G Mat ns (J(C)) if 
and only if p(\A\) < 1. 



5. A generalization: Positive semirings 

A semiring S with zero is positive if it is partially ordered by a relation 
=<! such that is the least element and for all x,y, z £ S the inequality 
x =<! y implies that x * z =<! y * z and z* x =<! z*y (see, e.g., [27]). Any 
idempotent semiring is positive with respect to the canonical partial 
order. The semiring of nonnegative real numbers M+ with the usual 
addition and multiplication and the ordering ^ provides an example of 
a nonidempotent positive semiring. 

Most of constructions and results of this paper hold for positive 
semirings with some minor changes. 

The semiring of square matrices Mat nn (5) over a positive semiring S 
is positive with respect to the following ordering: A = (a^) =^ B = (bij) 
in M&t nn (S) iff aij =4 hj for all 1 ^ i,j ^ n. By definition, a partial 
unary closure operation * in a positive semiring with unity 1 satisfies 
the conditions x* = 1 © x x* = 1 © x* x (in particular, 0* = 1) 
and x* =4 U* for all x,y £ S such that x y, provided that x* and y* 
are defined. In R + , for example, x* = (1 — x)" 1 if ^ x < 1 and x* is 
undefined otherwise. 

Weak and strong interval extensions of a positive semiring are de- 
fined similarly to the idempotent case. 

Suppose S is a positive semiring, I(S) is its weak interval extension, 
A € Mat nn (i"(5)), B G Mat ns (i"(5)), and the closure of the matrix 
A = [A, A] is A* = [A*, A*]. Consider the set S(A, B) of all solutions 
to the equation X = AX © B such that X = A*B, where A G A, 
B G B, and the algebraic solution A*B of the same equation. Then 
the following theorem holds: 

THEOREM 6. The interval matrix A*B G Mat ns (/(S")), considered as 
an element of I(M&t ns (S)), contains the set S(A,B) and is the least 
interval of I(M&t ns (S)) with this property. In particular, bounds of the 
interval A*B belong to E(A,B). 
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The proof is similar to that of theorem 4. Note that this theorem 
implies both theorem 4 and the theorem of W. Barth and E. Nuding 
cited above. 
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